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Abstract 

We study the expectation value of the energy momentum tensor during thin shell collapse for a 
massive, real, scalar field theory. At tree-level, we find thermal, Hawking-type, behavior for the 
energy flux. Using the Schwinger-Keldysh technique, we calculate two-loop corrections to the tree- 
level correlation functions and show that they exhibit secular growth, suggesting the breakdown of 
the perturbation theory. 


2 


I. INTRODUCTION 


Black holes are one of the few laboratories for studying quantum gravity. Indeed, a commonly 
viewed challenge for any candidate theory of quantum gravity is to provide an understanding of 
the quantum nature of black holes, and in particular their radiation. 

Hawking’s semi-classical calculation [l| shows that black holes radiate with a thermal spectrum. 
This has lead to black holes being viewed as objects to which thermodynamic quantities can be 
associated with, such as temperature and entropy, which satisfy relations analogous to the standard 
laws of thermodynamics. ^ The thermodynamic perspective on black holes is incredibly attractive 
and has inspired many questions regarding, for example, a microscopic understanding of its entropy, 
or whether it implies a holographic interpretation of gravity [^. However, the thermal spectrum 
of black hole radiation also poses a puzzle in the form of the information paradox j^. Our goal, 
in this paper, is to revisit black hole radiation and study quantum loop corrections to Hawking 
radiation that are often dismissed as negligible. 

Hawking radiation [l| (see also j^) is a quantum effect that can be seen in a gaussian (non¬ 
self-interacting) quantum field theory on a collapse background. ^ In this paper, we begin by 
reproducing Hawking radiation in a slightly different setting. 

The background that we would eventually like to consider is a star that is static due to an 
internal pressure until some moment of time, after which the pressure is switched off and collapse 
begins. To our knowledge the best model of such a situation is pressure-free collapse of a spherical 
perfect fluid—Oppenheimer-Snyder collapse [^. For simplicity reasons, however, instead of a ball 
of the perfect fluid, we study Hawking radiation in the case of a massive thin shell collapse. We 
will study particle creation during Oppenheimer-Snyder collapse elsewhere. 

The background in question is described in detail in the section |TT1 It has three phases, see figure 
m The phase labeled I in the figure is when the shell is kept fixed at some radius R{t) = Rq by 
an additional force. Phase H is a highly tuned stage of collapse if one wishes to respect spherical 
symmetry. In fact, if Rq is large enough, because of tidal forces any perturbation violating the 
spherical symmetry will grow in this phase. Phase HI, however, describes the final stage of the 
collapse, where, without rotation, spherical symmetry is restored (if it was lost in the second 
stage) because all multipole momenta are radiated away 10| , ll[ . This feature of the last stage of 
collapse is because of the peculiar properties of the horizon—the surface of infinite red shift. This 
phenomenon is also at the core of the “No Hair Theorem”. Hence, this stage is universal and does 
not depend on what is happening during phase H. 

It should not be surprising that a non-stationary gravitational background field creates particles. 
The production rate is then generically non-stationary and non-universal. But at the final stage of 
the collapse we have a stationary and universal free fall of the shell. Hence, it is natural to expect 
a stationary and universal particle production rate. Then the question reduces to the following 


^ While these laws were in fact suggested before Hawking radiation was proposed (see reference 0] and references 
therein), they can only be made sense of if black holes radiate. 

^ For reviews see, for example, references i, a and Ql- 
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Figure 1: The collapse of a spherically symmetric thin shell as is seen by an outside Schwarzschild observer. 
We divide the collapse into three phases: I, II and III 


one: What is the spectrum of the particles produced at the final stage of collapse? In other words, 
we expect that the free Hamiltonian of a theory on such a time-dependent background should also 
depend on time. But it is obvious that the free Hamiltonian can be diagonalized before the start 
of collapse. Its ground state defines the initial conditions for the problem that we study here. 
Moreover, we expect that it is also possible to diagonalize the free Hamiltonian at the final stage 
of collapse due to its stationarity. But one can guess that the harmonic functions which do such 
a diagonalization before collapse do not coincide with those performing the diagonalization at the 
final stage of collapse. The technical details supporting this general discussion are presented in the 
section imi 

In the section IIVI we calculate the time evolution of the expectation value of the stress-energy 
tensor. We find that before collapse has started the energy flux is zero. But at the final stage of 
the collapse it is given by a thermal stationary energy flux. We perform the calculation beyond 
the geometric optic approximation in four dimensional space-time for a massive real scalar field 
and include all spherical harmonics. 

The above discussion is true for a gaussian theory. In section El we study what happens if self¬ 
interactions, in particular a 4 >‘^ interaction, is turned on. The common wisdom is that quantum loop 
corrections should not change the picture described above in a substantial way. This is probably 
true for UV corrections. However, in section El we show that pertnrbative IR corrections grow 
with time and, if one considers a long enough period of time they can even dominate over the tree- 
level contribution. This is, in fact, a well known phenomenon in non-stationary condensed matter 
theory 0, 0: UV modes essentially do not feel the presence of the background field, but the 
behavior of IR modes reveals a change in the background state of the quantum field theory under 
consideration. In the concluding section we explain how these corrections modify the Hawking 
flux. 

As far as we are aware this is the first study of the effect of interactions on black hole radiation. 
For other work where various other corrections to the standard Hawk ing radiation picture are 
considered see, for 


example, 0, 0, 0, 0, 0, 0, 0, 0 and 0. 
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A. On the origin of secular growth in non stationary situation in perturbation theory 


In section |Vl we show that that the perturbative corrections to the Keldysh propagator of 
scalars on a collapsing shell background grow secularly. In this subsection we explain why this is 
not surprising and is to be generically expected in non-stationary situations. 

Suppose that one would like to find the time evolution of the expectation value of an operator 
O: 


(O),. (t) = (w 


ftQ Q Jig * ftg dt' H(t') 


4' 


( 1 ) 


The operator O could, for example, be the stress-energy tensor. Here H{t) = Ho(f) -|- V{t) is the 
full Hamiltonian of the theory; while T denotes time-ordering and T is anti-time-ordering; to is an 
initial time; IT) is an initial state. The initial value of the operator (O) (to) is known. We also 
assume some kind of covariant UV regularization. 

After transforming to the interaction picture, we get [1^ 


= /t 


S^t, to)S^oo, t) S{qo, t) T [Ooit) S{t, to)] 


T 


5^(cx), to) T [Ooit) S{oo, to)] 


T 


( 2 ) 


where S{t 2 ,ti) = Te * Ooit) and Vb(t) are the O and V operators in the interaction 

picture. The first step in equation ([2|) is an application of the Baker-Hausdorff formula. To perform 
the step on the second line of ([2|), we insert the unit operator 1 = 5^(oo,t) 5(oo,t), which allows 
us to extend the original time-evolution from to to t, 5(t, to), and back, 5^(t,to), to that which 
goes from to to future infinity, S(oo,to), and back, S'^(oo,to). Finally, we place the operator Ooit) 
on the part of the time contour that is going forward. 

Now let us slightly change the problem. Suppose we adiabatically turn on an interaction term, 
V, after time to- Then the state |T) does not change before to- In this case, we can rewrite the 
expectation value (l2|) as follows: 


it) = (T|5i(oo,-oo)r[Oo(i)5i„(oo,-oo)] 


T 


(3) 


An important question is if one can take to to past infinity, to —oo. If the answer is positive, then 
the remaining problem is if this can be done within perturbation theory or non-perturbatively. A 
well-know situation where this can be done within perturbation theory is when the free Hamilto¬ 
nian, Ho, does not depend on time and is bounded from below. Moreover, |T) should coincide with 
its true ground state, |uac), so that Ho \vac) = 0. In the calculation that follows, we also assume 
that the interaction term is adiabatically switched off at future infinity—after the time t- 

In fact, if the state |uac) is the true vacuum state of the free theory, then, by adiabat¬ 
ically turning on and then switching off the interactions, we cannot excite such a state, i.e. 
(^vac\S^ ioo,—oo)\ excited state) = 0, while |(uac |S'^(oo, —oo)| uac)| = 1. Thus, the dependence 











5 


on tQ disappears at every order of perturbation theory and we arrive at expressions that only con¬ 
tain time-ordering: we obtain an expression that can be calculated with the use of the standard 
Feynman diagrammatic method: 

(o>.. (*) = E ( vac |5^(oo, —oo)| state^ {state \T [Oo{t) S{oo, —oo)]| vac) 

state 

= ^uac |s^(oo, —oo) vac^ (uac |T [Oo(^) 5'(oo, —oo)]| uac) 


{vac\T [Oo{t) S{oo, —oo)]| vac) 


(4) 


(uac |5'(oo, —oo)| vac) 

In the first line in equation (j3|), we have inserted the unit operator 1 = 'Ylstate \ state) (state|, where 
the sum runs over a complete basis of states. In the second line, we have used the fact that |uac) 
is the only state in the sum that gives a non-zero contribution. 

If, however, |'I') is not a vacuum state, or Hq is time dependent (and, hence, cannot be di¬ 
agonalized once and forever), or Hq is not bounded from below, the above technique cannot be 
used and we have to directly evaluate the expression given on the right handside of equation Q. 
Moreover, in principle, the dependence on to can remain and it frequently does remain within per¬ 
turbation theory and sometimes even nonperturbatively—after resummation of leading corrections 
from all loops. The well known situation when the dependence on tg does disappear after such 
a resummation is when the system reaches thermal equilibrium from a nonstationary initial state 
[l^ . [l^ . This happens in flat spacetime without background fields, when Jig is time-independent 
and is bounded from below. A good question, however, is whether there is a stationary state in 
the presence of an eternally acting background field—such as for example global de Sitter space, 
constant electric held or eternal Schwarzschild black hole. Let us explain how the dependence on 
tg reveals itself. 

The standard method to calculate an expectation value such as that given in equation ([3]) 
is the so-called Schwinger-Keldysh diagrammatic tec hniq ue, where both S' and S^ have to be 
perturbatively expanded under the quantum average [1^, [1^. In such a case every vertex in 
the diagrams carries either a “-I-” or ” index, depending on whether it comes from S or S^. 
Furthermore, after Wick contractions every held, e.g. (/>, is described by a matrix of propagators, 
whose entries are given by 


n_+(Xi,X2) = (0 (Xi),^(A2)), B^_(Xi,X2) = (,^(A2)(/)(Ai)), 

B++(Xi,X2) = (T<)>(Xi)</>(X2)), B._(Xi,X2) = (T(/)(Ai)</-(X2)). 


(5) 


Propagators with time or anti-time-ordering or without orderings appear depending on whether 
the held 0 comes from S or S^. From their dehnitions, it is clear that these propagators obey the 
relation B^_-|- B_= B^^ -|- B _ 

After Keldysh rotation (l^ . [l^ . one of the entries of the propagator matrix is set to zero. 
The other entries are given by the retarded, advanced and Keldysh propagators. The retarded 
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propagator is given by 


{Xi,X2) = e (Atu) [D.+ iXi,X2) - D+_ iXi,X2)] 


— 0 {Ati2) [(l){Xi), (j){X2)\ , Ati2 —ti— t2- 


(6) 


While the advanced propagator, is conjugate to the retarded one. Since the commutator [•, •] 
is a c-number, at tree-level these two propagators do not depend on the state—they define the 
spectrum of excitations in the theory. 

The last entry of the propagator matrix is the Keldysh propagator: 

D^iXi,X2) = ^[D.+ iX,,X2) + D+.{XuX2)] = ^ mX,),(f,{X2)}), (7) 

where {•, •} is the anti-commutator. The Keldysh propagator is sensitive to the time evolution of 
the background state. This is an important fact for finding the perturbative secular growth that 
we are after in this paper. 

To explain the origin of the secular growth, consider the simple example of a spatially homoge¬ 
neous but non-stationary situation. ^ Then, in the case of a real scalar field cf) we have the following 
harmonic expansion: 



a^gkit)e'^^ + c.c. 


where the harmonic functions, gk{t) solve the corresponding Klein-Gordon equation in a time- 
dependent spatially homogeneous background field in flat space-time. It is not hard to see that if in 
equation ([71) we take the average with respect to an arbitrary state respecting spatial homogeneity, 
then: 


/' 


(t^xe (ti,x; t2, 0) = 


2 + yi“* 


Skih) + Hkiti) Hkitk) + c.c., (8) 


where we have used the fact that in spatially homogeneous situations 


If in such an expression we average over the vacuum, then, = 0 = Note that 

all the standard tree-level calculations of particle fluxes (such as e.g. Hawking flux or Schwinger’s 
electric current) are done with the use of such Wightman correlation functions as dSj) with ^at = 

0 = (afca_4- 

Furthermore, at tree-level and remain constant, even if the average is done 

over an arbitrary state, because all the time dependence of the creation and annihilation operators 


^ Of course in the case of spherically symmetric collapse background the details will be different, but conceptually 
the phenomenon is the same. The details for the collapse situation are presented in the section |V| Here we describe 
a simpler sitnation to give a flavor of the physical meaning of secularly growing quantum corrections. It is also 
worth stressing that in standard condensed matter situations one considers a non-stationary initial state rather 
than a background field. 
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is absorbed into the harmonic functions gk{t)- Note that the appropriate quantity to define as 


particle density is 


at a; 


9l{t)9k{t), rather than {ataj:'^ 

the form of the correlation function in equation ([ 8 ]). While for plane waves there is no difference, 
for the exact harmonics in background fields the difference is obvious. The same comment also 
holds for the anomalous quantum average However, below, for simplicity, we refer to 

particle density and to {a^a_^ as the anomalous quantum average. We hope that 
this does not cause confusion. 

If one turns on interactions, the behavior of the Keldysh propagator, ([ 8 ]), changes drastically. 
In particular, and (^a^a_j:'^ are generated, even if they were zero to begin with. That of 

course is to be expected: the presence of a non-zero anomalous quantum average means that, in a 
strong field background, the ground state of the theory is changing from past to future infinity once 
one turns on self-interactions. The same happens with the particle number density. In principle, 
they will depend on ti and t 2 simultaneously. However, since we are only interested in the leading 
growing terms, in the limit t = {ti -\-t 2)/2 —)• oo and \ti — 12 \ = const, we can neglect the difference 
between ti and t 2 - In fact, in theory, the two-loop sunset diagram correction to the Keldysh 
propagator is given by ([ 8 ]) with ^ 

(t) oc // dts (it4 6 (^k + qi + q 2 + qs) gu {h) gl (* 4 ) [ d^Qj gqj {h) (U), 


-this can be seen, in particular, from 


k 


to 






(t) oc -A^ f dt3 f dt4 6 (^+qi + q2 + qs) gl (ts) gl iU) [ d^qj g^^ (ts) g^. (t4) 

■/£q ■/£q j = \'^ 


(9) 


if the initial state is chosen to be the ground state of the free Hamiltonian at past infinity: 

\ground) = 0. Here to is the moment after which A is adiabatically turned on. In spatially 
inhomogeneous situations the formulae are quite different from the present one, but conceptually 
the phenomenon is the same. Also if the initial state is not a ground state then ([2]) gets modified. 

These formulae are also valid in the stationary situation, when the free Hamiltonian is time 
independent and is bounded from below. But it is not hard to see that if gp{t) oc then 

after change of integration variables from and ^4 to T = (ta -|- 14)/2 and t = — t^, the integral 

over r leads to a 5-function establishing energy conservation 5 [uj{p) + uj{qi) + ^( 92 ) + ^(gs)]. This 
5-function appears in the limit t — to —>■ 00 . On mass-shell the arguments of the 5-functions, 
establishing the energy and momentum conservation, cannot be simultaneously zero. Hence, in 
the absence of a background field, i.e. in a stationary situation, energy-momentum conservation 
forbids the change of level population and of the ground state. Then the ground state \ground) is 
actually a true vacuum state, |uac), of the free theory in question. Which agrees with the discussion 
before equation Q. 

However, in the presence of a background field there is no energy conservation (or the energy 
is not bounded from below), because the system is not closed. In equation Q this fact reveals 


Note that the bubble diagram correction to the propagator, ~ A, does not introduce corrections that grow with 
time because the corresponding contribution is local, i.e. is not sensitive to the position of the external legs. 
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itself via the presence of the exact harmonics gp{t), which are drastically different from the simple 
exponential. Therefore, the expressions for the level population and for the anomalous average are 
not zero at loop order. 

Furthermore, in the presence of a background field the harmonics are functions of physical 
momenta. For example, in the constant electric field, E, background harmonics are functions of 
k + eEt and, hence, are invariant under the simultaneous compensating shifts of k and t (see e.g. 

231, 241 for the detailed discussion). In de Sitter space, the physical momentum is ke~^^, where 
H is the Hubble constant. Hence, the harmonics are invariant under compensating shift of time 
and rescaling of the momentum (see e.g. 2^ for the detailed discussion). Finally, in the collapse 
background at future infinity and in the vicinity of the horizon, harmonics are functions of oo , 
where uj is the frequency and rg is the corresponding Schwarzschild radius of the collapsing body, 
as we will see in section UlI HI Because of these symmetries, after the change of integration variables 
from ta and t 4 to T = (ts + t4)/2 and t = in equation Q, one can get rid of the dependence 

of the integrand on T for low enough external physical momenta. In such a case one finds that 
(t) (X {t — to) and (t) oc (t — Iq), where (t — to) comes from the integration 

over T. The coefficient of proportionality here can be interpreted as a part of the collision integral 
that is responsible for the change of the population numbers from the background field and its 
backreaction on the initial ground state. One important complication with respect to the standard 
condensed matter theory situations is the presence of the anomalous quantum average. In standard 
situations (see e.g. [l^ and [l^) the ground state of the theory does not change. In the latter cases 
the only quantity that receives growing corrections, if the initial state is chosen to be nonstationary, 
is the particle number density and the formulae are a bit different from those in ([9]). 

In the previous paragraph, we have assumed that the background field is always on—from the 
past to the future infinity. However, if the background field is turned on at some time t*, then we 
can put to before t*. In such a situation, {t) oc (t — t*) and (t) oc A^ (t — t*). 

In either case we have the secular growth of the loop corrections. Thus, after a long enough time 
of evolution we obtain A^ t ~ 1 and loop corrections become of the order of tree-level classical 
contributions. The corrections, thus, have to be resummed in all loops. ^ Such a resummation has 


been done in the case of the constant electric field 


Sitter space massive scalar field theory in [2^ , [2^ , 


jack grou nds in 231, 12^ a nd in the case of de 
3, H (see also 0, 13 and Q, H)- 


In the present paper we start the same kind of study in the case of black hole radiation. 


® Note that this is a signal of the breakdown of perturbations theory, which has a clear physical meaning. The 
resummation which has to be done is the standard procedure in condensed matter theory non-stationary situations 
[l^ . [l^ . Apart from summation of the leading loop IR contributions it allows us to define the correct time 
evolution of the particle population numbers. In our case the main technical and conceptual complication is the 
presence of the anomalous quantum average, which makes this example similar to de Sitter space quantum field 
theory [^ . 
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II. THE BACKGROUND GEOMETRY 


In this paper, we consider a spherically-symmetric, massive, thin shell that is kept fixed at 
radius r = Rq (by an additional force) until t = 0 when it is released and collapses in free-fall— 
see figure [TJ Because of spherical symmetry, by Birkhoff’s theorem, the geometry is given by the 
Schwarzschild metric outside the shell and flat space inside: 


ds^ = 


dtl - dr^ - r^dn^, 


r < R{t) 


(l - r > R{t) 


dO^ = dd^ + cos^ Odip'^, 


( 10 ) 


where R{t) is the radial coordinate of the shell, which before the start of collapse is R{t < 0) = Ro] 
rg/2 is the ADM mass of the shell and t (t_) is the time coordinate outside (inside) the shell. We 
assume that Rq > Vg and, before the start of collapse, the shell is close to its Schwarzschild radius, 


i.e. |i?o - rg\ < rg. 


Note that, ideally, we would like to consider the collapse of some compact stellar object such as 
say a neutron star, whose radius is much larger than its Schwarzschild radius. However, here we 
assume that \Ro — rg\<^rg. In fact, in general, it is not possible for a stellar object of radius smaller 
than 3rgj2 to be stable. We consider such a phenomenologically unrealistic simplification, however, 
to make analytic headway. In particular, such a simplification allows us to find the behavior of the 
harmonic functions on such a background. 

The background that we consider is described by two metrics m that have to be matched at 
the position of the shell 3J, [3^. Before collapse, this is easily done because the radius of the shell 
is time-independent, i.e. on the shell r = Rq and dr = 0, and we have that 


= 



t < 0. 


( 11 ) 


After the start of collapse, we assume that the shell is in the free fall, hence the metric on the shell 
is ds‘^ = — R?‘{t) dVi?. Comparing this to the metric inside the shell, 



2 



= 1 . 


Furthermore, comparing the metric on the shell to the metric outside the shell gives a relation 
between t and proper time of the shell r, 




( 12 ) 


where R = ^. Using our assumption that \R — rg\ <C Vg, noting that i? / 0 and ^ oo as 
t oo, we neglect the right handside in comparison to the left handside and integrate the above 
relation to find the shell’s trajectory from the point of view of the outside observer: 


R{t) 



Rq - rg 



(13) 









10 


From the point of view of the inside observer, given our assumption that \Rq — rg\ rg, the shell 
collapses at constant speed v [3^, where 


dR{t-) 

dt- 


R 

\/l + ^2 ’ 


and R{t-) 


Rq — ut-. 


(14) 


Now, identifying R{t) = R{t-) and using equations (fT3|) and ([HI, we find the relation between t_ 
and t: 




Rq - rn 


1 — e 


oo. 


(15) 


As t —)• oo, the shell approaches its Schwarzschild radius R{t) —>• rg. On the other hand, ~ 
{Rq — rg) /i/ < oo is the moment of internal time when the horizon is formed. Beyond this moment 
of time there is no relation between t and In section imi we use the relations between t and t- 
before the collapse and at its final stage, relations m and (|15|) respectively, to find the behavior 
of the free scalar harmonics in these regimes. 

In deriving the behavior of the harmonics during collapse, we also use another simplifying 
assumption that we will now explain. Gluing the two metrics in (|10p to each other with the use of 
the shell’s energy-momentum tensor, we obtain the relation 3^ (see also 3f)l|): 


^ = mVi + r^ - —, 

2 2R' 


(16) 


where M = const defines the energy-momentum tensor of the shell’s matter content. This rela¬ 
tion expresses the fact that the total energy of the shell, rg/2, is the sum of the kinetic energy, 
M\/ 1 -KR 2 (recall that r is proper time rather than coordinate time), and its potential energy, 
M‘^/2R. To avoid a bounce, when R = 0, i.e. to have the actual collapse, we must have that 
rg > 2M. However, we in fact make a stronger assumption, namely that rg 3> 2M. Neglecting the 
difference between R{t) and rg, from equation (I16p we obtain 


R: 




\2M 


R 


3> 1 and 1 / ~ 1, i.e. the shell is almost 


le point of view of the inside observer. While 


Hence our assumption that rg ;§> 2M implies that 
light-like at the final stage of the collapse even from t 
we use 1 / ~ 1, we keep v explicitly in most of the formulae presented below to indicate how the 
situation would be different without this assumption. 

In the region outside the shell, it is convenient to also use tortoise coordinates 


r* = r -I- r„ log-1 

Vg 

In terms of these coordinates the trajectory of the shell at the final stage of the collapse is 

R^{t) ^ Rl - t + {rg - Rq) (l - e 


(17) 


where Rq = Rq + rg log ~ 1) j i-®- the shell’s motion becomes light-like as t —)• oo—see figure 

El 


















11 



Figure 2: The collapse of a massive, thin shell in tortoise coordinates. We again divide the collapse into 
three phases: 1, 11 and 111 


III. FREE HARMONICS 


The theory that we study on the background of a collapsing shell is the real scalar held theory 
with (j)‘^ interaction: 


S= d‘^x^/\i\ 




(18) 


In this section, we discuss the free theory case A = 0. 
We expand the held cj) as follows: 


4> {t, r, 0 , ^ Yi^ri {0, (1)1 {t, r), 


Ln 


where (0, (p) = Yi^n (0, pc) are real spherical harmonics. Given this expansion and the back¬ 
ground given by equation (| 10 ll . the action takes the form 


+ E (2' + 1) / “t f - (l - i) 

Varying this action, we obtain the equations of motion, 

- 0“^ + 7n? + {r^i) = 0, r < R{t) 

dl - dl* -h (l - ^) [rn^ + -h (r^z) = 0, r > R{t) 
and the boundary conditions, 

(1)1 R{t) - 0 = (1)1 R{t) + 0 , 




2 \ i2 


■ (19) 


( 20 ) 


If dt\ 

dR 

^ , / dt \ ^ , 



dR 

( 1"a\ ^ , 

l(d 

dt 

dtfpi - 

r=R{t)-0 

_ £9 

r 

dt 

- (1 - j dr(l)l 


( 21 ) 


r=R{t )+0 
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The second equation above relates the normal derivative of the scalar field across the shell. As 
t —7> oo, the derivative normal to the shell outside it becomes the derivative with respect to the 
retarded coordinate u = t — r*. 

The equations of motion outside the shell can be understood in terms of confluent Heun equa¬ 


tions 37 h 39I| and there has been recent progress in understanding the properties of solutions to 
these equations 40!,l4]|. However, for our purposes it is sufficient to consider approximate solutions 
of these equations in various stationary regions of space-time. In this way we obtain estimates of 
the leading effects/contributions to Hawking radiation. 

The reason we can find approximate solutions to equations (I20p and m lies in the fact that 
in the stationary situation we can separate the variables t and r. Furthermore, the potential. 


U{r) = 


2 l{l +1 ) r, 


r > Ro, 


( 22 ) 


in the second equation of (|20l) is vanishing as r ^ Vg and approaches a constant, m?, as r ^ oo. 
Therefore, in the vicinity of the shell—note the relation \R{t) — r^l <C Vg —the harmonics can be 
approximated by plane waves. Meanwhile, at asymptotic spatial infinity, r ^ oo, we have that 
r ~ r* and the harmonics are again plane waves. The relation between the harmonics in the 
vicinity of the shell and at spatial infinity can be found by solving the scattering problem across 
the potential barrier (1221) . In the next subsection, we will use the above arguments to find the 
approximate harmonics. 

Using the equations of motion, the free Hamiltonian of the massive scalar theory can be rewritten 
as 


Ho{t) = ^(2Z + I) 

i 


Jo sm(6») 


{dtcPif 


1 


4>idt 


(23) 


In the collapsing shell background, using the metric given in equation (lIOp . this Hamiltonian is 


Ho{t) = E 

l^n 


R(t) 



0 






{Mif 


4'idt4)i 


+ E 

l,n 



dr 4^1 


r=R{t)—Q 



drCpl 


r=R{t)+0 


(24) 


The last term in this expression is the contribution of the field values on the shell itself. This 
Hamiltonian defines translations along the time coordinate t. The corresponding Cauchy surfaces 
are depicted in the figure [3l 


A. In-harmonics before collapse 

During the first stage, when the shell is stationary and is yet to collapse, we can find harmonics 
that diagonalize the free hamiltonian, hence providing a sensible definition of particle number. The 
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singularity 



Figure 3: Penrose diagram of the collapsing shell background. The curved, thin lines depict Cauchy surfaces 
with respect to the Schwarzschild time t. The grey region represents the shell interior. 


state with respect to which the flux is found is defined in terms of these harmonics, which we call 
in-harmonics. 

As is usual in canonical quantization, we expand the scalar field in terms of a basis of harmonics 



auj,i,ndth^i{r,t) + h.c. 


(25) 


where, at this stage, u simply labels the harmonics and “h.c.” stands for hermitian conjugate. The 
harmonics provide a basis for solutions of the Klein-Gordon equation (I20p . Note that we 

only expand in harmonics with oj > m because only these modes are oscillatory at spatial infinity. 
We explain this point further below equation (I3UI1 . The harmonic functions that we wish to work 
with are such that they diagonalize the Hamiltonian when the shell is stationary. We refer to the 
corresponding ground state, which is annihilated by all annihilation operators, = 0, as 

the in-state. 

The canonical commutation relations for the field cj) (x, t) are 

[<P{x, t),TT{y, t)] = {x - y) (26) 


and all other commutations relations vanish. Assuming, 

~ 2705;;/ dnn' ~ ^ 


(27) 


we obtain the following condition on the modes from the canonical commutation relations 

foo u, , r_ _ 1 

/ 75 - h^,iit,r)dthAiit,r') -h.c. = i <5^^^ (x - x') . 

1 ^ Jm L ’ J 


( 28 ) 
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We stress that x are the spatial coordinates. 

In a stationary situation we can separate the dependence of hu),i{t,r) on t and r. Then during 
the stage before collapse, the positive energy states can be represented as 


K,i{t,r) = h^^l{r)e 


hoj,i{r) 




where from equation (fTTTl it follows that co- = uj /Given the decomposition above, we 
can identify the label ui in equation ([2^ with the energy of the mode hi^^t(t,r) at past infinity. 
Inside the shell, by solving the first equation in (12011 . we find that 





for r < Rq and t < 0. 


(29) 


Note that we only require Bessel functions of the first kind, which lead to finite harmonics at 
r = 0. ® From equation (j29ll . we can see that a; is bounded from below: uj- > m or io > m_ = 
m — Tg/RQ. If the mass of the scalar is small, Au) can be approximated by y/lr —see appendix 
s Note that it is not necessary to make this assumption and we could simply keep A^] in the 
expressions without solving the normalization condition (j28l) . However, we will use Aui ~ for 
convenience. 

Similarly, from the second equation in (j2nh . we find that, before t = 0, the harmonics outside 
the shell are given by 


J A^e \r - i?o| < rg, 

where k = \/uA — m?. 

Now we can see that when uj > m, the harmonics oscillate as r* —>■ oo. Hence, in such a 
case we have the continuous spectrum. But when m_ < u < m the harmonics decay or grow 
exponentially as r* —>■ oo. After dropping the exponentially growing mode, we are left with modes 
that correspond to a discrete spectrum that oscillate between r = 0 and the turning point of the 
potential (I22p . Prom the gluing conditions across the potential (|22n one can estimate the allowed 
values of the discrete energy levels Ui. We leave the study of these discrete modes for the future. 
For our present purposes, it is important to know that there are only finitely many such states in 
the well, which has a depth of order m? rg/R^. Hence, these states play a minor role in the effects 
that we consider here. In particular, since they do not correspond to running states at spatial 
infinity, they do not contribute to Hawking radiation at tree-level. Moreover, if we take m?rg/RQ 
to be small enough—we are in fact working with small m —there will be no discrete states at all. 

The absolute value of the undetermined coefficients in equation (1301) can be found using the 
normalization condition (I28p . However, since we are only interested in the behavior of the har¬ 
monics in the vicinity of the shell, we will concentrate on coefficients and B^. The effects that 



6 


Note that these Bessel functions lead to the diagonalization of the free Hamiltonian in empty flat space-time in 
spherical coordinates. 
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we are looking for appear in the vicinity of the collapsing shell. As seen from spatial infinity, or 
elsewhere, these effects simply receive grey body factors due to the potential barrier (1221) . 

We find coefficients and using the boundary conditions (f2T]) . During the first stage, 
when R{t) = Rq, the boundary conditions reduce to 


h, 


(jJ^l 


Rq — 0 


= h. 


'LJ,l 


Rq + 0 




In appendix [HI we show that the above equations imply that 

„•/+! / r \ 


a^ = b: = ^ 

v2a; 


(-a 




J +0 


1 - 


Rc 


3/4N 


(31) 


(32) 


Let us consider the free Hamiltonian for these modes (see reference 4^ for a similar discussion). 
In a generic situation, the free Hamiltonian, expressed via harmonic functions, takes the following 
form; 


Ho{t) = 


Ln 


diddio' 


£ui,uj'A A h.C. 


(33) 


where 

^0. 


J. 


= [ 

Jo 

r 

hio'At) = 

Jo 




1 


dr 


sin(0) 

VW\ 

sin(0) 


9^^dtK i{r,t)dtK>^i{rR) - —hl i{r,t)dt ^A9\ ,i{r,t) 

visl L 

9 ^^dtK,i{r,t)dtK',i{r,t) - ^h^^i{r,t)dt \^/\g\ g^^dthu',i{r,t) 

VlffI L 


Thus, generically the free Hamiltonian is not diagonal because of the presence of terms with 
u ^ ijj' and especially because of the presence of a non-zero However, given suitably chosen 

normalizable modes, in stationary situations £ and J become constants because of the factorization 
of the dependence of the harmonic functions on t and r: h^^i{r,t) = /i^^;(r). In such a case 

one can find harmonic functions, satisfying 

V\i\ 


d?’ h^',i{r) = 0, 

Sine? 


f 


dr 


5 “ K',i{f) = -S{uj - uj'), (34) 

sm u ’ UJ 


which lead to the vanishing J' terms. Note that the normalization above is chosen so as to be 
consistent with condition (|28p in stationary situations. 

The stage before collapse, considered in this section, is indeed a stationary situation, but we 
have discussed the behavior of the modes for certain ranges of r. To be able to find the free 

Hamiltonian, we integrate over r and hence, in principle, we need the modes for all r. However, 
here it is sufficient to assume that there are normalizable modes that satisfy equation (I34p — these 
modes will necessarily have the form given in equations (1291) and (I30p for the appropriate regions. 
With the normalization of the modes given in equation (1341) . it is simple to verify that the free 
Hamiltonian is 

duj 


r 

Hoit <0) = ^ 

l,n 


27r 


-UJ 






(35) 


and is indeed diagonal before the start of the collapse. But during collapse, the in-harmonics no 
longer diagonalize the Hamiltonian. This is a sign that there is particle creation. 
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Figure 4: An illustration of the Cauchy problem for harmonics after collapse. 


B. In-harmonics during the late-stage of collapse 


In this section, we find the behavior of in-harmonics as t — > oo. The change in the behavior 
of the in-harmonics from past to future infinity due to the non-stationarity of the background is 
what causes particle creation. 

We assume that z/ ~ 1 and neglect the difference between u and 1. The reason for this is 
as follows: Given a harmonic function inside the shell, see figure 01 we would like to find the 
harmonic function outside the shell by solving the gluing condition (j21|l . To find the form of the 
harmonic function behind the shell one has to solve the Cauchy problem with initial values given 
at t = 0 = t_. The initial condition is given by the values of the harmonic functions behind the 
shell and just outside it and their derivatives at t = 0 = f_. Then, at a point inside the shell 
which is causally disconnected from the region outside of the shell at t = 0, e.g. a point labeled 
“o” in figure (j4]), the value of the harmonic is simply given by equation (I29j) . However, in order to 
use the boundary condition to find the harmonic functions outside the shell, we need to solve the 
Cauchy problem for the region inside the shell, grey region in figure 01 that is within the domain of 
influence of points outside the shell on the initial Cauchy hypersurface. This appears to be a rather 
complicated problem. However, if we take ~ 1 then the domain of dependence of points inside 
the shell are always inside the shell, hence the grey region can be neglected and the matching can 
be done straightforwardly, as we will see below. 

Assuming z^ ~ 1, the harmonics inside the shell are given by equation (I29p . where uj- = 
u / Y^l — Tg/RQ. Since we are only interested in the modes near the exterior of the shell, assuming 
\R{t) — r^l <C Tg, we neglect the potential ([2^ . The solution is again as in the case before collapse 
and can be given in terms of Fourier modes. However, as the matching at the shell is more 
complicated with the solution in this form, we instead keep the solution general. In all we have 
the following solution, as t —>■ +oo. 




1 


+5'd^,z(u). 


r 


r < R{t), 

r > R{t), \r - R{t)\ < Vg, 


( 36 ) 
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where u = t — and v = t + are outgoing and ingoing null coordinates, respectively. Since the 
shell is light-like at the late-stage of collapse, the u-dependent part of the harmonics outside the 
shell is not affected by the collapse. This can be easily understood by noting that the boundary 
conditions are continuity of the harmonics and their derivatives normal to shell. Hence, guj,i{v) is 
unaffected by the collapse and, from equation ([30]) . we have 

(37) 

We now find the u-dependent part of the harmonic outside the shell by imposing the boundary 
conditions (|2ip at the shell. Note that we are assuming that ~ 1, which means that the shell 
collapses at approximately the speed of light even from the point of view of observers inside the 
shell. In terms of the picture in figure HI this means that we are assuming that the grey region 
is infinitesimally thin. In this approximation the harmonics inside the shell are given by Bessel 
functions, as in the static case, but close to the shell they also have a constant contribution that 
is equal to the value of the u-dependent part of the harmonics outside the shell. Therefore, taking 
this constant piece into account, for the positive frequency modes the continuity of across the 
shell gives 

■ (38) 

For brevity we denote R{t) as R, but the dependence of i? on t should be remembered. Using 
equation (fT7|) . as t —> oo the outgoing null coordinates evaluated at the shell are approximately 


[u]r=R ~ - {Rl + rg- Ro). 

Hence, the continuity relation reduces to 

fuj,i [2t - {Rq + rg- i?o)] ~ ~ ® 

Therefore, using relation (1151) . we find that 


—iuj-t— 


(39) 


(40) 


, , I u + Rf^+rg-Rg 

-lui- ^ ° 1-e 


r hui,i{r, t) K, R{u)Jij^^ - m? R{u) 

where g^j^iiy) is given in equation ((371) and 

R{u) = rg 

As t —)• oo, the boundary condition for the derivative of the modes at the shell, m, becomes 


— Tn u+Rq+Vq-Rq 
^ y ^ 2rn 


+ guj,i{v), (41) 


(42) 


Since, as t —> oo, R —> Vg, we have used ^ « 1 to simplify the right-handside of the above 
equation. In AppendixjCl we show that this relation is satisfied by (1411) if we neglect the difference 
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between v and 1. In the case where u is much less than 1, our solution for the u-dependent part of 
the harmonic will be modified and will be more complicated. 

It is worth stressing at this point that the behavior of the harmonic function very far from the 
gravitating center (at spatial infinity) is the same both when the shell is static and when the shell 
is collapsing. That can be expected on general physical grounds because of locality. Hence, while 
the harmonics have dependence on null coordinates near the shell surface, this will not be the case 
at spatial infinity. This is expected because these modes are harmonics of a massive scalar field 
which propagates to future time-like infinity and not future null infinity. 


IV. HAWKING RADIATION (TREE-LEVEL) 


We will, first, reproduce the standard thermal radiation using the harmonics constructed in the 
previous section. In particular, we calculate the energy flux at the final stage of collapse. It is 
technically easier and physically more appropriate to calculate it just outside the shell, r > Vg, and 
then continue across the potential barrier (1221) . In fact, the flux is created just outside the black 
hole. Thus, we find the flux as t —>■ oo in the vicinity of the shell, |r — rg\ <C Vg. 

In the proximity of the shell, the energy flux is given by 


J 


{r^rg,t) = J sin 6 d6 d(p t (r, t) ps — J sin 6 d9 dip {r,t) 


'S2 ' ' JS2 

Using the expression for the energy-momentum tensor of a scalar field 


J{r-rg,t) = ^(2/ + I) 


U / ? 


(44) 


(45) 


where 


= r 


j 


= r 


m 

oo 


doj 

2^ 


Jm 


duK^i (r, t) duK,i {r, t) + c.c. 

dvKArit) dji^^i{r,t)+C.C. 


(46) 


In equation (I44p . the normal ordering is understood in the standard way in the presence of the 
background fields. In particular, we understand (: Ty :) as a subtraction from (Ty) of the same 
expression but calculated for the case of a shell that is eternally static at r = i?o- In this case, 
since the harmonic functions in the vicinity of a static shell are given by Fourier modes in u and 
V, from equation (|45p . the flux vanishes. 

In order to find the flux due to the collapsing shell we use the harmonics outside the shell at 
the late-stage of collapse given in equation (l4T]l . 




uj-rg 


u — uq 

gi uirg e ‘■‘’'3 



r 


a 





'/2uj 


—iujv+iLj[RQ—Ro{l—rg/Ro) 

^ 5 


( 47 ) 
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where uq = rglog{Ro/rg) — Vg. In order to obtain this expression from equation (jH]), we use 
R{u) ~ Tg, z/ ~ 1; neglect the difference between r and Vg in the denominator and neglect m? in 
comparison with = cj^/ (1 — Vg/Ro), because Rq ps Vg and a; > m for the continuous part of the 
spectrum. We also use the limiting form of the Bessel function of the half-integer index with large 
argument. 

In order to calculate the flux, we now simply need to insert the harmonic given in equation (j47|] 
into equation (I45jl and evaluate the integrals. However, in order to make contact with the original 
derivation of Hawking radiation [l| and for simplicity, it is appropriate to re-expand the expression 
on the right handside of equation (|47l) in terms of the modes given in equation (I30p . Because the 
u-dependent part is unchanged we only have to re-expand the u-dependent part: 




2 

— cos 

OJ 


Tr{l + 1) 


- uj-r„ 


^lUJTge 


u — uo 
' ‘2ra 


doj' 




2Tr^/^J 


; a{uj^uj ) e 


/\ —iuju 


(48) 


where = a(w, |a;'|) and = a{LJ, — |a:'|) are proportional to the seminal Bogoliubov coef¬ 
ficients for Hawking radiation. 

The implicit expression for a{u},u}') given in equation (I48p is in fact only valid for f > 0. There 
is also a contribution to a{uj,uj') from before collapse. However, these contributions will be of the 
form S(cj — oj'). Since we are only interested in the effects due to collapse, we will ignore these 
terms in a{uj,uj') and, from equation (Hsp . find 

/ 

J u* 


a 




2(1-^ 

Ro 




cos 


LJ 


7r{l + 1) 


- UJ-Vn 


due^' 


'9 ruj'u 
^ 1 


(49) 


where u* = —Rq- Because of the rapid oscillations of the integrand in the lower limit of the u 
integration compared with the upper limit, the exact position of n* is not very relevant and we can 
extend the range of integration to the real line. Moreover, as we have eluded to before, since we are 
interested in the late-stage flux generated by the collapsing shell, we have ignored contributions 
to a (w,a;') from stages I and H of the collapse (see figure [2]). Evaluating the u integration on the 
real line, it is straightforward to show that 


a 


Rq 




cos 


u 


tt{1 + 1) 


- oj-r. 


Using the above results, the expressions for the and ^ are 


gTTO; Vg ^2iu'rg log {u^Vg ) p ^ , 

(50) 




1 - 


Tg \ 2 r°° duj 

Ro) Jm 27r 


UJ 


and 


jd) 

On, 


1 

f dio' 

r duo” 

uo'uo” 

Jm 27r J 

CN 

J\uj''\>m 2 vr 

■\/\uo' uj”\ 


aiuj^uj ) a (uj^uj ) ^ ^ • 


(51) 


(52) 


(0 


The integral over uj in the expression for Ju can first be evaluated as follows: 

/■oo J . . i 


TOO d(log (a;rg)) 2i(uj'-uj")r„ log (corn) Zco', 

J log (mrg) 2 -k i \ ) 
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where, since oj-Vg i:|> 1, we have replaced cos^ 


7r (l+l) 


- OJ-ra 


by 1/2 in the integral—other con¬ 


tributions lead to terms that decay as powers of 1/u in equation (j52p . The integral above is over 
the real half-line, hence 


[ ^ a{ui,ui') a* 2rg |a;'| |r(2ia;Vg)|^ 5(0;'— w")-|-regular term, 

Jra V *^0 / 

1 

/ \ 2 

(54) 


= 2 TT (1 —^ I n{—uj') 6 fuj' — Lo") + regular term, 
' ^0/ 


where 


n(uj) = 


sign(uj) 

g47rrgaj _ ^ ' 


(55) 


The “regular term” on the right handside of equation (1541) is of the form p.v. j , hence for 

large u its contribution to Ju \ ([52]) . is negligible. 

Substituting equation ([541) into Ju'^ and using n(—w) = n{u) + sign{uj), we obtain the following 
expression for the flux: 



(56) 


Therefore, the total flux is 


J{r^rg,t) = ^{2l + 1) 


v(2/+i) 

RoJ ^ Jm 27r 


Lo n{u}). 


(57) 


This is the black body radiation appearing in the vicinity of the collapsing shell. The flux at 
infinity is modihed due to the mass of the field under consideration and is multiplied by the grey 
body factor, due to the potential barrier (1221) . 


J {r ^ oo, t) 



(58) 


where at large I the grey body factor ^ can be estimated to behave as ^( 2 l-i)!! example 


i+i 


reference 


m- 


If we naively set m = 0 in equations ()48]) - (|M]l . we encounter divergent integrals at w = 0. 
However, in the massless case one has to use 


U-UQ 


hui,i{r,t) K. y/W^Jij^i{uj-rg)e H-e 




(59) 


instead of equation ()17)l . It is not hard to show that because of the peculiar behavior of the Bessel 
functions at small values of their argument, (x) ~ all oj integrals in equations (|48]) ~ ([58]) 

are convergent at their lower bounds. Then the calculation of the flux continues in the same way 
as for the massive case and we obtain equation (|58p where m is set to zero. 
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V. LOOP CORRECTIONS AND THE SECULAR GROWTH 

In the previous section we used the two point Wightman function to calculate the expectation 
value of the stress-energy tensor. As motivated in the introduction section, here we show that loop 
corrections to this Wightman function grow with time. 

The interaction term for theory in the collapsing shell background has the following form 

oo 

^ ^ (^) / (60) 

R{t) 0 

In the limit t ^ oo, the second term is exponentially suppressed because oc Hence, we 

can neglect this term. Furthermore, in this limit, we can neglect the difference between R{t) and 
Vg as the integrals are regular at r = r^. 

As mentioned in the section [U the one-loop, bubble diagram correction to the Keldysh propa¬ 
gator does not introduce a secularly growing contribution, hence we consider the next loop order. 
The two-loop sunset diagram contribution to the Keldysh propagator can be expressed as follows: 

Uo^+2(1,2) = I 

h,m\,l2,m2 

\j^ijj\,li,ni\u>2,12,^2 T ^hh ^mim2 <^(^1 ^ 2 )] (H Tl) ^aj2,i2 (^2; ^ 2 ) 

,ni\u]2,h,n2 (^) h{ti,ri) 

hu)2,l2 (t2U2) + ll.C.j, (61) 

where we have used the definition of the Keldysh propagator in terms of U_|_and D _(_ propagators, 

equation d?]). Also, from equation (j5|). 

poo J 

D+-{l,2) = J2 ■:^K,iiti,ri)h^^i{t2,r2)Ypn{0i,(pi)Ypn{02,(P2), (62) 

1,771 ^ 

roo j 

D-+il,2) = J2 ■7^h,^Ah,ri)hl, i{t2,r2)Ypn{0i,(pi)Ypni02,V’2) ■ (63) 

l,m ^ 

Neglecting the difference between ti and t 2 and letting t = + ^ 2 ) 72 , it can be shown that the 

coefficients N and K in the Keldysh propagator take a simple form, 

f f - - 

Nu,,i,n\io',i\n'{t) = ^ // dtsdti {r3r4A drsdr4h^^i{r3,t3)hl,,i,{r4,t4) 

^ doj ■ 

Y{l,n,l' j (64) 
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and 


^ t t3 ^ 

Kaj,i,n\uj',i',n'it) = - Y j ^^3 J dU JJ (rsrif dr^ dr^ {/i* ^(rs, ta)t4) + (w, / o w',/')} 


^0 ^0 


/*CO . 

F(/,n,/',n') fj / (65) 

j^l^m 


where 


Y{l,n,l',n') = ^ {I,n,li,ni,l 2 ,n 2 ,l 3 ,n 3 ) {I',n',li,ni,l2,n2,k,n3) , 

h,2,3',ni,2,3 

{I, n, Ii,ni,l2,n2,l3, ^ 3 ) = J dQ Yi „^(9, cp) Yi^ ^i {d,p) Yi^^n 2 (d-,p) ^( 3,^3 {e,p). 

Note that we use the real basis of spherical harmonics. 

It is worth noting that N and K from (1641) themselves do not receive UV divergent 
contributions—this can be seen from simple power counting even in flat space-time. However, 
the Keldysh propagator has the standard UV divergence due to integrations over oj. Our main 
interest, here, is the contribution from IR modes to N and K. Hence, we assume some suitable 
UV renormalization and assume that all coupling constants acquire their physical (renormalized) 
values. 

We start with the consideration of V^j,z,n|i.j',z',n'(^)- Our goal is to single out the largest contri¬ 
bution to N and K in the infinite future limit. This limit corresponds to uj e u)' e 0 

for fixed oj and oj '—we will keep only the leading terms in this limit. Note that there are no 
large contributions to N and K during the stationary stage before collapse because of energy 
conservation—this was explained in the introduction section. Hence, we can put to before t = 0— 
before the onset of collapse. Furthermore, as we argued at the end of section m very far from the 
collapsing shell the behavior of the harmonic functions practically does not change after the start 
of the collapse. Hence, there are no growing contributions from the region of large r. Therefore, 
the harmonic functions in equation (fMl) can be approximated by (1471) . i.e. ps h^^i{u) + huj,ii'v), 
where huj,i{v) is the n-dependent part of huj,i{r,t) given in equation (H71) . while huj,i{u) is its u- 
dependent part. Note that we still keep 00 as the upper limits of integration over r 3 ^ 4 , because the 
integrals are rapidly converging—in fact the integrals converge much faster than if the modes at 
spatial infinity were used. Moreover, the largest contribution to V)^,/,n|a;',/',n'(i) comes from the re¬ 
gion of integration over t 3 and t 4 where t 3 S> Vg log {vg uj) and t 4 S> rg log (r^ u'). In this region we 
can neglect the dependence of h^^i{u 3 ) and /i*, ;/(u 4 ) on U 3 and U 4 , respectively. The n-dependent 
parts of hi^^i{r 3 ,t^) and h*, ;,(r 4 ,t 4 ) lead to subleading contributions. 
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Using the above simplifications, we find that 


^u),l,n\u]',V,n' (^) 


2A2 


3r2 


OJUJ' 


Rc 


1 ^9 \ ^ 

1 — COS 


TT (/ + 1) 


- W_r„ 


cos 


TT {I' + 1) , 

--- ijj_r, 


-'9 


rt pt poo poo 

Y{l,n,l',n') / dts / dt^ / r|dr 3 / r|dr 4 

Jrglog(rga;) JVg log{rg uj') Jr„ Jr„ 


3 

n 

i=i' 


duj , 


2 ^ huj,,igir3,t3)hl,^i.{r4,U). ( 66 ) 


We can again express the harmonics as sums of u and z;-dependent parts, hence 

roo . _ _ POO , p_ 

J ir3,t3) K. ;. (r4, t4 ) « J {U3 ) hi. ^iYuA) + h^. {V3 ) hi. ^i. {v^) 

+ Kg,lY'^3)Kg,lg{VA) +K.^iYV3)K^ • (67) 

Since, hu)^i{u) is proportional to a rapidly oscillating cosine function, the last two terms in the 
above expression are negligible. The first contribution, can again be expanded in Fourier modes, 
as in the previous section. Using equations (|151) and (|M|) . 


dujj , 


in, t3)hl.^i. (r4, t4) ~ { 1 - ^ ) 4 


dbji 


RoJ 


{ n(-ujj) 


\ -iujj(u3-ui) 


+n{ijjj) +e + subleading terms, (68) 

where n(a;) is defined in equation (l55|) . 

We can now substitute equation (1681) into the expression for N^^i^n\uj',i',n'i't), (t66|) . Changing the 
integration variables from t 3 and t 4 to T = (t 3 + 14)/2 and t = t 3 — t 4 , we obtain 
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(69) 


where Ar = r 3 — r 4 . Note that the lower limit of integration over dT is irrelevant for the leading 
term that grows with time. We have also extended the limits of integration over r to ±oo because 
the rapid oscillations of the integrand makes their exact position irrelevant as t —>■ oo. 

Now we can see that the integrand of the dT integral does not depend on T itself. Hence, the 
leading contribution is A ~ t. As explained in the introduction, it is also possible to see how in 
a stationary situation this contribution would vanish. In such a case, all the exponentials would 
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have the same sign in their exponent, hence the r integration would give (5(a;i +UJ 2 + W 3 ), which is 
zero given that iOj are positive. 

Let us now show that the same growth appears in K. Making the same approximations in 
equation ()65p as is made for the calculation of N, we arrive at the following expression: 
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Then, using equation () 68 p and performing the same change of integration variables from t 3 ,t 4 to 
T and r as above, we get that 
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We again see that the integrand of the integral over T does not depend on T itself. Hence, here 
the leading contribution is also K XH. 

For the massless case the situation is similar, but instead of using the harmonic of the form 
given in equation (1471) . we must use h given in equation (l5^ to evaluate 

/ OO _ 

■^K,l{^3,t3)h^^l{u,t4). (70) 

However, as before J;_,_i(w) is negligible in the small uj region, so the only difference with the 
massive case is that the uj integrals run from 0 to oo. Therefore, we obtain the same type of 
behavior of N and K for massless fields. 


A. Loop corrections to the Hawking radiation 


Above we have shown that loop corrections to the two-point functions are not suppressed in 
comparison with the tree-level contribution. In fact, even if A^ <C 1 after a long enough time of 
collapse A^ t ~ 1. It means that perturbation theory breaks down. In particular, higher loops 
bring higher powers of }?t to N and K. Thus, one has to do the resummation of all leading loop 
corrections. This kind of resummation was done for the case of a strong electric field and massive 
scalar field theory on de Sitter space-time in 23], 24], 25], 26], 27], 28], 29], 32], 33|. For the 


present case the resummation will be done elsewhere. 
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Of course it is too early to draw any conclusions before this resummation has been done. 
However, based on the intuition gained from previous examples, we make a few remarks on the 
physical consequence of these growing corrections. Their presence means that after resummation 
one will get some specific time-dependence Ku),i,n\u)',v which does not 

necessarily have to be linear. Furthermore, on general physical grounds we expect the presence 
of a stationary final state for the case under discussion. However, we expect that due to the loop 
nature of the corrections in question this stationary state will be reached at much latter times than 
the stationarity of the Hawking tree-level flux is reached. 

In any case the presence of N and K means the modification of the Hawking flux. In fact, in 
the case when iF = 0 and N ^ 0 the equation for the flux is as follows: 
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(71) 


It is not hard to also restore the contribution of K. Now because of the loop corrections we have 
that 


^1^2,1,— 6{u}i ^^ 2 ) + 2 ( 1 ' 2 ) 

The above tree-level contribution is obtained when = 0. Note that due to the orthogonality of 
the spherical harmonics only diagonal in I and n components of contribute to the 

flux. Thus, one can see that the presence of N does lead to a modification of the Hawking flux. 


VI. CONCLUSIONS 


In this paper, we show that, in the case of Xcj)'^ theory, loop corrections to the Hawking flux are 
not suppressed given a sufficiently long time—in other words, perturbation theory breaks down. 
In order to make a definitive conclusion about the fate of Hawking radiation and quantum black 
holes it is necessary to consider the resummation of the leading corrections from all loops, as done 
for the case of a strong electric field and massive scalar field theory on de Sitter space-time 23(], 
24l |. 25|, 2^, 29(1 (see 28(] for a review). We leave the resummation of the corrections, for 

the present case, for the future. 

The presence of unsuppressed loop corrections leads to a modification of the Hawking flux. 
Furthermore, while the tree-level flux {T[) is universal, loop corrections to it, i.e. N and K, are 
not universal. Their final value after resummation and as t 00 will, in general, depend on initial 
conditions. At the same time, in some approximation they completely characterize the state of the 
quantum field theory in the vicinity of the black hole and affect the spectrum of the radiation. In 
any case, in this paper we have demonstrated that it is inappropriate and premature to draw any 
conclusions about black hole evaporation and puzzles associated with it, such as the information 
paradox j^, based on a study of the tree-level contributions to the energy-momentum tensor. 
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Appendix A: Normalization of the harmonic functions 

In this appendix, we determine the normalization of the harmonics inside the shell before collapse 
using the condition from the canonical commutation relation. As we remark in section IHH it is not 
necessary to fix the normalization in order to find onr results—we simply fix the normalization for 
convenience. For the modes inside the shell, given by equation (j29|) . the left handside of condition 
(128|) reduces to 

lIs (' ■ %) ') ■''+4 ’■') 

2 (pr) {pr') , 

where we have made use of the following substitntion for uj: 

p = — m?. 

When mrg<c(l — ,we can use the normalization identity for Bessel functions whereupon 

we obtain ^ 

(Al) 

7rr^ 

Hence, 

Aoj ~ v^. (A2) 



^ We deduce that must be p-independent in order to obtain the form of the normalization identity for Bessel 

functions. 
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Appendix B: The boundary condition for the harmonic functions on the shell before collapse 


We fix the coefficients and of the harmonic modes in the vicinity of the shell, ()30p . by 
imposing the boundary condition (I21|) on the shell. First, consider the continuity of the modes at 
the shell, (IfiTT) . From jRo — Vgj ^ Vg, we deduce that 
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Hence, up to lower order terms in ^1 — , the continuity of the (pi across the shell gives 

T b - ^ ii, - h W , 
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which, in particular, implies that B^j = . Furthermore, the second boundary condition in (j.SlIl 

for the radial derivative of pi, again up to lower order terms in gives 
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Appendix C: The boundary condition for the harmonic functions on the shell during the 

late-stage of collapse 

In this appendix, we show that the solution for the harmonic function outside and in the vicinity 
the shell, m, which was found in section IIIIBI by imposing continuity at the shell also satisfies 
the boundary condition on the normal derivative, equation (143 jl . 

In section [III Bl we added a constant term to inside the shell because of the approximation 
1 / ~ 1. In the boundary condition for the normal derivatives of the harmonics this term simply 
cancels the u-dependent term in outside the shell. This makes sense because the boundary 
condition is a condition on derivatives normal to the shell, namely the u-dependence of the har¬ 
monics outside the shell. In what follows, we will concentrate on the non-constant piece inside the 
shell and the u-dependent term in the harmonics outside the shell. 















28 


Using the harmonic function inside the shell, given in equation (1361) . the left-handside of equation 
(H3l) is 
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where we have used the expression for t_, CSl), as t —)■ 00 . The right-handside of equation ()^3]l is 
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where we have used equation 
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and the behavior of u on the shell as t —> 00 , ([39]). Using the above relation and equation ([39]) . 
the right-handside of equation (|13]) simplifies to 
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Comparing the above equation to equation (ICll) . we conclude that the boundary condition on the 
normal derivatives is satisfied \i v = \. 
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